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摘 要 : 而 砂岩 区 是 黄土 高 原水 土 流失 最 严重 的 地 区 之 一 ,其 植被 生长 状况 对 该 区 控 表 
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土壤 侵蚀 和 维护 生态 平衡 


起 着 重要 作用 。 本 文 基 于 CASA 模 型 和 Rclimdex 1.0 分 别 计算 了 2001 一 2021 EIEI E KIERRE I NPP 和 


18 个 极端 气候 指数 ,采用 


趋势 分 析 、 相 关 分 析 、 随 机 森林 重要 性 排序 . 残 差 分 析 等 研究 1 


磊 砂岩 区 NPP 时 空 变化 及 


其 对 气候 因素 的 响应 特征 ,并 量化 了 气候 因素 和 人 类 活动 对 而 砂岩 区 NPP 的 相对 贡献 。 结 果 表 明 :(1) 2001—2021 
年 础 砂岩 区 NPP 变 化 呈 显 著 增 加 趋势 ,但 未 来 82.5% 的 区 域 NPP 将 由 增加 趋势 变 为 减少 趋势 。(2) 年 际 尺 度 上 flt 
砂岩 区 NPP 与 年 均 气 温 .年 总 降水 量 .极端 强 降水 指数 主要 呈正 相关 ,与 冷 夜 日 数 TNI0P .气温 日 较 差 DTR 呈 负 相 
关 。 季 节 尺 度 上 ,春季 均 温 以 及 暖 夜 日 数 增 加 有 利于 NPP 增 加 , 且 存 在 滞后 影响 。 夏 季 暖 则 日 数 增加 不 利于 植被 


生长 , 且 NPP 对 夏季 暖 凡 日 数 存在 3 个 月 的 滞后 响应 。 夏 季 极 端 强 降 水 有 利于 NPP 增 加 ,而 夏季 
, 且 NPP 对 持续 干燥 日 数 CDD 存 在 3 个 月 的 滞后 响应 。(3) 气候 因素 和 人 类 活动 共同 促进 了 而 砂岩 区 NPP 增 加 ， 


类 活动 贡献 率 为 60.40% 。 


关键 词 : 础 砂岩 区 ;植被 净 初 级 生产 力 (NPP)， 极端 气 候 ; 


植被 净 初 级 生产 力 NPP 是 植物 光合 作用 所 固 
定 有 机 质 总 量 中 扣除 自 养 呼吸 后 的 剩余 部 分 , 即 绿 
色 植物 在 单位 面积 和 单位 时 间 内 积累 有 机 物 的 数 
量 "。NPP 作 为 植被 生态 系统 中 物质 循环 和 能 量 流 
动 的 基础 ,与 生态 系统 功能 紧密 相关 ,反映 出 植物 
群落 在 自然 环境 条 件 下 的 生产 力 , 是 评价 植被 生态 
系统 可 持续 发 展 的 重要 生态 指标 ,是 判定 生态 系统 
碳 源 / 汇 的 主要 因子 ”。 和 气候 因素 和 人 类 活动 是 影 
Iie] NPP 的 两 个 重要 驱动 因素 "。 已 有 研究 表明 ,人 
为 土地 利用 管理 是 全 球 植被 变 绿 的 重要 驱动 因素 ， 
贡献 了 三 分 之 一 的 植被 叶 面 积 增 量 ,尤其 在 中 国 和 
印度 **。 而 无 直接 土地 利用 变化 影响 的 地 区 ,各 生 
态 系统 植被 变化 在 大 多 数 情况 下 可 归 因 于 气候 系 
统 变量 的 趋势 "。 在 大 多 数 区 域 植被 同时 受到 气候 
因素 和 人 类 活动 的 影响 ,影响 程度 因 地 区 而 异 。 因 
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量 不 利于 植被 生 


露 区 和 覆 沙 区 气候 贡献 占 主导 地 位 ,气候 因素 相对 贡献 为 62.13% 和 60.06% ,覆土 区 以 人 类 活动 贡献 为 主导 ,人 


人 类 活动 ; 残 差分 析 


此 ,量化 气候 变化 和 人 类 活动 对 植被 变化 的 相对 贡 
献 , 有 利于 区 域 植被 的 进一步 保护 .建设 和 管理 , 同 
时 对 理解 区 域 植被 动态 变化 成 因 ,制定 合理 的 区 域 
生态 恢复 策略 具有 重要 意义 。 

气候 因素 中 温度 和 降水 是 影响 植物 生长 最 直 
接 和 关键 的 气候 因子 ,许多 研究 分 析 了 平均 气候 
因素 对 NPP 变 化 的 影响 1, 一般 来 说 ,在 干旱 和 半 
干旱 地 区 ,降水 增多 对 水 分 限制 区 植被 产生 积极 
影响 。 在 北半球 高 纬度 等 热量 限制 区 ,温度 升 高 
减轻 了 温度 对 植被 生长 的 限制 ,进而 有 利于 NPP 
增加 后 和 。 与 平均 气候 相 比 ,极端 气候 对 NPP 的 影 
响 较 少 被 关注 ,但 是 NPP 对 极端 气候 的 响应 可 能 
为 敏感 '"。2000 一 2009 年 干旱 导致 全 球 陆 地 净 初 
级 生产 减少 "1, 暴 雨 和 长 时 间 降 雨 也 会 限制 植物 呼 
吸 和 光合 作用 六 。 另 一 方面 ,人 类 活动 也 是 影响 植 
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被 NPP 变化 的 重要 原因 ,生态 恢复 成 效 在 不 同 区 域 
可 能 存在 差异 。 在 我 国 黄 土 高 原 ” ,横断 山区 "等 
区 域 人 类 活动 是 NPP 增 加 的 主要 原因 ,表明 生态 修 
复工 程 对 改善 植被 NPP 起 着 重要 的 作用 。 三 峡 库 
区 实施 的 生态 修复 工程 有 效 缓解 了 气候 变化 的 不 
利 影响 ,主导 了 森林 恢复 ,而 中 下 游 地 区 快速 的 城 
市 化 则 导致 了 森林 退化 " 1。 虽然, 已 有 许多 研究 量 
化 了 气候 变化 与 人 类 活动 对 植被 NPP 变化 的 相对 
贡献 ,但 大 多 数 NPP 预测 模型 是 基于 降水 和 温度 等 
平均 气候 因子 ,同时 考虑 极端 气候 作为 气候 因素 的 
研究 较 少 。 因 此 ,本文 考虑 极端 气候 指数 以 量化 气 
候 因 素 和 人 类 活动 对 植被 NPP 变 化 的 相对 贡献 ,可 
在 一 定 程 度 上 减少 气候 因素 对 植被 变化 解释 的 不 
准确 性 。 

而 砂岩 区 位 于 晋 陕 ` 蒙 三 省 交界 处 ,是 黄河 中 
游 剧烈 侵蚀 区 ,其 基 岩 成 岩 程 度 低 .胶结 程度 差 ,无 
水 时 坚硬 如 租 石 , 遇 水 则 溃散 如 泥 ,该 区 域 水 土 流 
失 严 重 , 是 黄河 中 游 粗 泥 沙 的 主要 来 源 区 ,生态 环 
境 极为 脆弱 "55。 目 前 ,针对 看 砂岩 区 NPP 的 研究 较 
少 , 平 均 气 候 、 极 端 气候 及 人 类 活动 对 该 区 NPP 变 
化 的 影响 尚 不 明晰 。 因 此 ,本 文采 用 CASA 模型 佑 
FAW A DX NPP , oP Hr Al WD X. NPP 时空 变 化 及 其 
对 平均 气候 和 极端 气候 的 响应 ,并 结合 随机 森林 重 
要 性 排序 结果 ,考虑 极端 气候 指数 进行 多 元 残 差分 
析 ,量化 气候 因素 和 人 类 活动 对 看 砂岩 区 NPP 变 化 
的 相对 贡献 ,以 期 为 而 砂岩 区 应 对 气候 变化 下 的 植 
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1983 


被 恢复 工作 提供 相应 的 科学 依据 ,为 而 砂岩 区 生态 
环境 建设 及 黄河 流域 高 质量 发 展 提供 参考 。 


1 数据 与 方法 


11 研究 区 概况 

tt RP 44 X fiz F 38° 42' ~40° 12' N, 108? 52' ~ 
111?35' E Z [B] ,属于 黄土 高 原 北部 , 陕 . 晋 、 蒙 三 省 
交界 的 鄂尔多斯 高 原 ( 图 1) ,面积 约 1.67x10' km’, 
地 势 西 北 高 东南 低 ,地 表 沟 蜜 纵横 ,土壤 侵蚀 严重 ， 
生态 环境 脆弱 ,是 黄土 高 原水 土 流失 最 为 严重 的 区 
IZ, MIARE PEE FY PON PR Be KBE E 
DX LUD DC ,根据 侵蚀 程度 又 分 为 剧烈 侵蚀 裸露 区 
和 强度 侵蚀 裸露 区 。 而 砂岩 区 年 均 气 温 6~9 % ,年 
均 降 水 量 315~442 mm ,为 温带 干旱 半 干 旱 大 陆 季 风 
气候 ,气候 干燥 ,降水 多 集中 于 夏季 , 占 全 年 降雨 量 
83 70% An 1. 2020 AF Ak RP zr DC TE HR TR RR ig TR 
05 61.796 ,林地 占 4.81% ,耕地 占 18.5996 ,其 他 土地 
和 覆盖 类 型 占 14.9%( 图 1) , EP zt CERCASI P] CR 
1.0 数据 来 源 及 处 理 

NDVI 数据 来 源 于 GEE 平台 MOD13Q1 数据 
(https://code.earthengine.google.com/) ,空间 分 辨 率 为 
250 m, 时 间 范 围 为 2001 一 2021 年 ,时 间 分 辩 率 为 
16 d。DEM 数 据 下 载 于 中 国 科学 院 地 理 空间 数据 云 
ASTER GDEM 数据 (https://www.gscloud.cn/) ,空间 分 
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图 1 而 砂岩 


区 2020 年 了 


上 地 利用 类 型 及 位 置 示意 图 


Fig. 1 Land cover types in 2020 and map of geographical location in the Pisha sandstone area 
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辩 率 为 30 m。 土 地 覆盖 类 型 数据 来 源 于 中 国 科 学 1.3.2 Sen+M-K 趋 势 分 析 本 文采 用 Theil-Sen Me- 


院 地 理科 学 与 资源 研究 所 (https:Wwww.resdc.cn/) o 
气象 站 点 数据 来 自 中 国 气象 数据 网 (https:/data. 
cma.cn/) ,为 了 保证 空间 插值 结果 更 加 精确 ,本 人 研究 
选取 了 而 砂岩 区 及 其 周边 17 个 气象 站 点 ,获得 
2001—2021 年 的 气温 、 降 水 量 和 日 照 时 数 数据 ,由 
于 太阳 辐射 站 点 较 少 ,采用 日 照 时 数 和 纬度 的 经 验 
公式 ”计算 大 阳 辐 射 。 极 端 气候 指数 从 ETCCDI 
(Expert Team for Climate Change Detection and Indi- 
ces) 所 提供 的 极端 气候 指标 体系 中 选择 了 10 个 极 
端 气温 指数 和 8 个 极端 降水 指数 ( 表 1) ,并 利用 
Relimdex 1.0 计 算 18 个 极端 气候 指数 。 最 后 统一 利 
用 ANUSPLIN 软件 进行 空间 插值 ,得 到 气象 栅 格 
1.3 研究 方法 

1.3.1 NPP 的 估算 本 文采 用 朱文 泉 等 "改进 的 
CASA EALAH TE IER DX NPP. CASA 模型 通过 植 
PO UR ACER SG-AA ARCA (AP AR) AG AE RU HIS (e) fih 
算 NPP ,计算 参数 少 ,并 已 得 到 广泛 应 用 ”” ,具体 
计算 方法 见 文献 ”2”。 


dian 非 参数 趋势 计算 方法 , RIC SERIE RP AK 
NPP 年 变化 趋势 。 该 方法 已 广泛 用 于 分 析 长 时 间 序 
列 数 据 的 趋势 变化 站 。 计 算 公 式 为 : 
NPP. - NPP, T 
IA e , 2001 <i<j<2021 
(1) 
式 中 ;median 表 示 中 位 数 函 数 ;6B 表示 NPP 的 变化 趋 
势 ; NPP, 和 NPP, 表示 不 同年 份 的 NPP 数 据 ; t, I t, 
为 时 间 序 列 。 当 B>0 时 ,表示 NPP 呈 上 升 趋势 ;68<0 
时 , 则 表示 存在 下 降 趋 势 。 并 采用 Mann-Kendall 检 
验 ( 简 称 M-K 检 验 ) 对 趋势 分 析 结 果 进 行 显著 性 
检验 。 
1.3.3 Hurst 指数 ”Hurst 指 数 用 于 预测 时 间 序 列 数 
据 的 未 来 变化 趋势 ,Hurst 值 ( 力 用 于 判断 未 来 变化 
趋势 ,0<H<0.5, 表明 未 来 变化 趋势 与 过 去 相反 。 
0.5<HE1, 表 明 未 来 变化 趋势 与 过 去 一 致 。 大 0.5， 
表明 时 间 序 列 为 随机 序列 ,未 来 趋势 不 确定 。 将 
NPP 变 化 趋势 与 Hurst 值 亚 加 ,以 确定 而 砂岩 区 NPP 
未 来 趋势 。 
1.3.4 相关 分 析 采用 Pearson 方 法 分 析 而 砂岩 区 
2001—2021 ^F. NPP 与 气候 因子 的 相关 性 ,计算 公 


B=median 


表 1 极端 气候 指数 


Tab.1 Definition of extreme climate indices 


ETT 代码 名 称 定义 单位 
极端 气温 。 TNIOP 冷 夜 日 数 日 最 低 气 温 < 10% 分 位 值 的 日 数 d 
冷 极 值 TX1op 冷 层 日 数 日 最 高 气温 <10% 分 位 值 的 日 数 d 

TNn 日 最 低 气温 极 小 值 。 ” 每 月 内 日 最 低 气温 的 最 小 值 iG 
TXn 日 最 高 气温 极 小 值 。 ” 每 月 内 日 最 高 气温 的 最 小 值 is 
极端 气温 。 TN90P 爱 夜 日 数 日 最 低 气 温 > 90% 分 位 值 的 日 数 d 
BRE 。 Txgop BAK 日 最 高 气温 > 90% 分 位 值 的 日 数 d 
TNx 日 最 低 气温 极 大 值 。 ”每 月 内 日 最 低 气温 的 最 大 值 c 
TXx 日 最 高 气温 极 大 值 。 每 月 内 日 最 高 气温 的 最 大 值 c 
其 他 气温 DTR 气温 日 较 差 最 高 气温 与 最 低 气温 的 差 什 c 
"OX GSL 生长 季 长 度 日 平均 气温 第 一 次 连续 6d 以 上 大 于 5 %C 至 日 平均 气温 第 一 次 (6 月 1 日 d 
后 ) 连续 6 d 小 于 5 %C 的 日 数 
极端 降水 ”PRCPTOT 雨 日 降水 总 量 十 日 (日 降水 量 =1 mm) 降水 总 量 mm 
指数 RX1day 1d 最 大 降水 量 每 月 最 大 1 d 降 水 量 mm 
RX5day 5 d 最 大 降水 量 每 月 连续 5 d 最 大 降水 量 mm 
R95P 强 降 水 量 每 年 日 降水 量 >95% 分 位 值 的 总 降水 量 mm 
ROOP 极 强 降 水 量 每 年 日 降水 量 > 99% 分 位 值 的 总 降水 量 mm 
SDII 年 均 雨 日 降水 强度 降水 量 =1l mm 的 总 量 与 总 日 数 的 比值 mm: d” 
CDD 持续 干燥 日 数 日 降水 量 < 1 mm 的 最 长 连续 日 数 d 
CWD 持续 湿润 日 数 日 降水 量 =1 mm 的 最 长 连续 日 数 d 


1238] ERES: URL RAAK DT BAY E EB DUAE ^ 7 BY 


Bloat Bon 

式 中 :R 为 x.y 两 变量 的 相关 系数 ;x 为 第 i 年 NPP 
值 ;y; 为 第 i 年 所 分 析 的 气候 指数 值 ;x 表示 多 年 
NPP 平 均值 ; y 表示 多 年 所 分 析 的 气候 指数 平均 
值 。n 为 样本 数 , 即 21。 

1.3.5 随机 森林 重要 性 排序 随机 森林 (RF) 是 一 
种 基于 分 类 和 回归 的 机 器 学 习 算 法 ,能够 用 来 分 析 
自 变 量 与 因 变 量 之 间 复 杂 的 非 线 性 关系 下。 本 研 
究 使 用 随机 森林 来 识别 各 个 极端 气候 对 而 砂岩 区 
NPP 的 影响 程度 ,利用 R4.2.1 中 的 RandomForest 包 
执行 RF 算法 ,得 到 增长 的 错误 率 平 方 均值 %Inc- 
MSE, B} Increased in Mean Squared Error (26) ,表示 某 
— AE RIS EROS EL POLI VER BEY S5 , 4o Inc MSE 
值 越 大 重要 性 越 高 ,用 以 直观 反映 各 极端 气候 指数 
对 目标 变量 (NPP) 的 重要 程度 。 

13.6 多 元 线性 回归 残 差 及 相对 贡献 分 析 多 元 线 
性 回归 残 差分 析 可 用 于 区 分 气候 因素 和 人 类 活动 
对 NPP 变化 的 影响 ,已 得 到 较为 广泛 的 使 用 。 
一 般 来 说 ,多 元 回归 残 差 分 析 选 取 平均 气候 变量 如 
总 降水 量 和 平均 气温 作为 气候 因素 ”” ,忽略 了 极 
端 气候 对 植被 的 影响 ,可 能 导致 对 气候 贡献 估计 的 
不 准确 性 ”。 因 此 ,本 研究 将 极端 气候 指数 考虑 进 
来 ,以 年 均 气温 \ 年 总 降水 量 和 极端 气候 指标 为 自 
变量 ,NPP 为 因 变 量 进 行 多 元 回归 残 差分 析 ,并 结合 
相对 贡献 分 析 方 法 ” ,计算 气候 因素 和 人 类 活动 对 
而 砂岩 区 NPP 变 化 的 相对 贡献 。 步 又 如 下 :(1) 以 
NPP 为 因 谈 量 , 年 均 温度 TEM ,年 总 降水 量 PRE AK 


1985 


端 气候 指数 incidex 为 自 变 量 ,建立 多 元 线性 回归 模 
型 (公式 3), 其 中 ,incidex 为 随机 森林 重要 性 排序 对 
fil» e X NPP 影响 程 度 最 高 的 极端 气候 指数 ,a,b， 
c,E 为 模型 参数 。 NPP. 为 逐年 NPP 预 测 值 ,表示 气 
候 因 素 对 植被 NPP 的 影响 。(2) 计算 NPP, (CASA 
模型 模拟 的 年 均 NPP 值 ) 与 NPP., 之 间 的 差 值 (公式 
4), Bl NPP, ,表示 人 类 活动 对 植被 NPP 的 影响 。 
NPP, =aXPRE+bXTEM+cXincidext+e (3) 
NPP,,, = NPP.,, - NPP. (4) 
(3) 计算 2001 一 2021 年 NPP_ 和 NPP， 的 线性 趋势 
3 , EH Slope( NPP... ) fll Slope( NPP, ), LA Be NPP 的 趋 
势 率 Slope( NPP，), 并 根据 表 2 计 算 气 候 因 素 和 人 
类 活动 的 相对 贡献 率 。 


2 结果 与 分 析 


2.1 NPP 模 拟 和 气象 插值 结果 评估 

分 别 在 CASA 模型 估算 的 年 NPP 和 MODIS NPP 
数据 集中 随机 抽取 100 个 像 元 , 共 获 得 2100(21 x 
100) 个 样本 ,表现 出 较 好 的 拟 合 结果 (图 2a) ,基于 
CASA 的 NPP 值 与 MOD17A3 NPP 呈 显著 相关 性 ,RR 
为 0.7068(P<0.01), 上 共有 一 定 的 统计 学 意义 。 同 时 ， 
ASSO ALIAS A EX NPP 年 均值 为 197.68 g Coma, 
TAS TR 5 xU TETE AE OU dE X EAR BE BS GERE 
岩 区 NPP 佑 算 范 围 一 致 。 整 体 上 ,本 文 基于 改进 的 
CASA 模型 估算 而 砂岩 区 植被 净 初 级 生产 力 精 度 

采用 交叉 检验 进行 气象 数据 插值 精度 误差 检 
验 , 选 取 而 砂 岩 区 空间 分 布 均匀 的 4 个 气象 站 点 作 
为 检验 样 点 ,剩余 的 作为 训练 样 点 进行 空间 插值 。 


表 2 气候 因素 和 人 类 活动 相对 贡献 率 计算 方法 


Tab.2 Methods for calculating relative contribution rates of climate factors and human activities 


lope(NPP,,.) 驱动 因素 驱动 因素 的 划分 标准 了 驱动 因素 的 贡献 率 /% 
slopel s 驱动 因素 
i ' slope(NPP..) slope(NPP,,,) 气候 变化 人 类 活动 
>0 CC&HA >0 >0 slope(NPP..) slope(NPP,,,) 
slope(NPP.,,.) slope(NPP.,.) 
CC >0 <0 100 0 
HA <0 >0 0 100 
<0 CC&HA <0 <0 slope(NPP...) slope(NPP 1s) 
slope(NPP.,,.) slope(NPP.,.) 
CC «0 20 100 0 
HA >0 <0 0 100 
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本 人 研究 选取 东 胜 、 保 德清 水 河 、 河 曲 4 个 站 点 作为 
检验 样 点 。 结 果 如 图 2 所 示 , 年 均 温 度 (R=0.77， 
P<0.01) 和 年 降水 量 (R=0.69 , P < 0.01) 与 实测 值 拟 
合 效 果 良 好 ,气象 数据 插值 精度 较 高 。 
2.2 而 砂岩 区 NPP 时 空 变化 

如 图 3 所 示 ,在 时 间 变 化 趋势 上 ,2001 一 2021 年 
间 而 砂岩 区 年 均 NPP 增 长 速率 为 5.78.a (P< 
0.001) ,NPP 总 体 呈 波动 上 升 趋势 , 且 而 砂岩 各 区 域 
NPP 波 动情 况 相 似 , 均 呈 显 著 上 升 趋势 。 其 中 ,覆土 
区 NPP 增 长 最 快 , 增 速达 7.02.a ,裸露 区 NPP 增 长 
速率 最 小 ,为 3.72.a 7 。 在 空间 分 布 上 ,2001 一 2021 
年 NPP 平 均值 与 覆土 程度 密切 相关 , IURI RPG 
NPP HEER IX (157.45 gC.m a) < EY, (203.57 
g Cem®- a") < H EX (217.38 g C* m?*a"), 

在 空间 变化 趋势 上 (图 4a), 2001—2021 4E 24 
96.94% 的 区 域 NPP 呈 增加 趋势 ,84.60% 的 区 域 呈 显 
著 增 加 趋势 。 不 显著 增加 的 区 域 (12.34% ) 主要 位 


于 裸露 区 。 结 果 表 明 ,2001 一 2021 年 间 而 砂岩 区 
NPP 明 显 增 加 ,植被 生长 状况 明显 改善 。 进 一 步 计 
算 Hurst 指数 , 75 E Alt WD A X. 2001—2021 年 植被 
Hurst (ELI 3-35 0.45 ,表明 未 来 NPP 的 变化 主要 
表现 为 反 可 持续 性 。 在 未 来 空间 变化 趋势 上 (图 
4b),82.50% 的 区 域 NPP 由 增加 变 为 减少 ,1.96% 的 
区 域 持续 减少 ,14.44% 的 区 域 持续 增加 ,1.1% 的 区 
域 由 减少 变 为 增加 。 这 一 结果 表明 ,植被 生长 未 来 
因 和 气候 变化 和 人 类 活动 易 发 生 退 化 ,为 了 促进 夏 砂 
岩 区 植被 长 期 稳定 生长 以 及 区 域 生 态 环境 可 持续 
发 展 ,巩固 当前 生态 恢复 工程 已 有 成 效 , 仍 需要 加 
强 对 而 砂岩 区 植被 的 保护 力度 ,实施 适合 该 区 域 的 
可 持续 性 发 展 的 植被 生态 恢复 工程 。 

2.3 而 砂岩 区 NPP 对 气候 因素 的 年 际 响应 

2.3.1 年 NPP 与 气温 的 相关 性 BIS ih, this 
区 年 NPP 与 气温 因子 的 空间 平均 相关 系数 绝对 值 
均 在 0.3 以 下 , 呈 弱 相关 关系 。 具 体 而 言 ,与 NPP 相 
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图 2 净 初 级 生产 力 与 气象 数据 精度 验证 


Fig.2 Accuracy verification of net primary productivity and meteorological data 
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图 3 2001 一 2021 年 NPP 时 间 变 化 趋势 及 空间 分 布 
Fig. 3 Temporal trends and spatial distribution of NPP from 2001 to 2021 
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(a) 2001 一 2021 年 变化 趋势 
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1987 


(b) 未 来 变化 趋势 
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图 4 2001 一 2021 年 NPP 空 间 变 化 趋势 及 NPP 未 来 变化 趋势 
Fig.4 Spatial trends of NPP from 2001 to 2021 and future trends of NPP 
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Fig. 5 Spatial distribution of correlation coefficients between NPP and extreme climate index in the Pisha sandstone area 


关 性 通过 显著 性 检验 面积 占 比 较 大 的 指标 为 冷 夜 
日 数 TN10P 气温 日 较 差 DTR 年 均 温度 TEM 与 生 
长 季 长 度 GSL。TEM、GSL 与 NPP 以 正 相关 关系 为 
主 ,表明 年 均 温 升 高 与 生长 期 延长 有 利于 看 砂岩 区 
植被 生长 。TN1I0P .DTR 与 NPP 以 负 相 关 关 系 为 主 ， 
表明 随 着 冷 夜 日 数 TN10P 的 增加 ,可 能 导致 植被 积 
温 不 足 ,影响 植 被 生长 。 而 气温 日 较 差 DTR 是 日 最 
高 气温 与 最 低 气 温 的 差 值 , 受 极端 气温 的 影响 , 当 
DTR 增 大 时 ,表明 昼夜 可 能 出 现 极 暖 或 极 冷 ,不 利 
于 植被 生长 。 其 中 ,TN10P DTR 与 NPP 呈 显著 负 相 


关 的 区 域 位 于 裸露 区 、 履 沙 区 以 及 部 分 覆土 区 南 
部 ,表明 相 较 于 覆土 区 ,极端 气温 对 裸露 区 和 有 覆 沙 
区 可 能 产生 更 消极 的 影响 。 此 外 ,上 暧 极 值 (TNx、 
TXx、TN90P、TX90P) 、 冷 极 值 (TNn、TXn ,TX10P) 5 
NPP 以 负 相 关 关 系 为 主 ,但 几乎 未 通过 显著 性 检 
验 。 总 体 来 说 , 极 暖和 极 冷 均 不 利于 而 砂岩 区 植被 
生长 。 

2.3.2 年 NPP 与 降水 的 相关 性 与 极端 气温 指数 相 
反 ,极端 降水 指数 与 NPP 表现 为 普遍 的 正 相关 (图 
6) ,上 且 相关 系数 较 高 ,通过 显著 性 检验 的 面积 更 
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图 6 MEDEX NPP 与 极端 降水 指数 相关 系数 的 空间 分 布 


Fig. 6 Spatial distribution of correlation coefficients between NPP and extreme precipitation index in the Pisha sandstone area 


大 。 年 总 降水 量 PRCPTOT .年 均 雨 日 强度 SDII 与 
NPP 呈 显著 正 相关 的 区 域 覆 盖 而 砂岩 区 绝 大 多 数 地 
区 。 不 同 强度 的 极端 降水 (RXlday、RX5day , R95P , 
R99P) 与 NPP 均 以 显著 正 相 关 为 主 , 且 空间 分 布 上 
基本 一 致 , 呈 显 著 正 相关 的 区 域 主要 分 布 在 裸露 
区 EU XANES. NPP 与 持续 湿润 日 数 
CWD 呈 显著 正 相关 的 区 域 主要 分 布 在 而 砂岩 区 秦 
土 区 和 和 覆 沙 区 的 南部 ,这 表明 短期 的 强 降水 有 助 于 
缓解 旱情 ,促进 裸露 区 植被 生长 ,而 持续 降水 可 能 
对 裸露 区 植被 产生 消极 影响 。 此 外 ,NPP 与 持续 干 
燥 日 数 CDD 平 均 相 关系 数 仅 为 0.04, 日 大 多 数 未 通 
过 显著 性 检验 。 总 体 来 说 ,人 研究 时 段 内 的 总 降水 量 
以 及 降水 强度 增加 对 全 域 NPP 均 产生 了 有 利 影 
We, CAE FH Tube XT T ETSI ,植被 生 
长 易 受 水 分 的 限制 ,降水 增加 促进 了 植被 生长 。 
2.4 而 砂岩 区 NPP 对 气候 因素 的 滞后 响应 

从 图 7 中 可 以 看 出 ,而 砂岩 区 NPP 主要 在 春季 


恒 日 数 TX90P 与 同期 NPP 表 现 为 显著 负 相 关 ,分 别 
存在 2 个 月 .3 个 月 的 清 后 效应 ,表明 夏季 上 暖 尽 日 数 
增加 不 利于 植被 生长 ,但 夏季 气温 日 较 差 减 小 有 利 
于 NPP 增 加 。 夏 季 RXlday、RX5day、CWD PRCP- 
TOT 均 与 同期 NPP 表 现 为 显著 正 相 关 , 量 存在 不 同 
程度 的 滞后 效应 ,表明 夏季 降水 增多 对 而 砂岩 区 植 
被 生长 可 产生 1~3 个 月 的 积极 影响 。 此 外 ,夏季 持 
续 干 燥 日 数 CDD 与 NPP 表 现 为 显著 负 相 关 , 且 存在 
3 个 月 的 清 后 效应 ,说 明 夏 季 干 旱 会 对 而 砂岩 区 植 
被 生长 产生 滞后 3 个 月 的 不 利 影响 。 
2.5 极端 气候 对 丰 砂 岩 区 年 NPP 变 化 影响 的 重要 
性 排序 

基于 随机 森林 重要 性 排序 测定 各 极端 气候 指 
BON AIL WD Ar X NPP 变化 影响 的 重要 性 ,结果 如 图 8 
所 示 , 气 温 日 较 差 DTR、 生 长 季 长 度 GSL、 雨 日 降水 
总 量 PRCPTOT 对 全 域 NPP 的 影响 在 气候 指标 中 排 
名 前 三 , 且 DTR 的 重要 性 远大 于 GSL、PRCPTOT。 


和 夏季 对 气温 和 降水 指标 产生 显著 的 响应 , 且 存 在 
一 定 的 滞后 性 。 春 季 , 冷 夜 日 数 TN10P 与 春季 同期 
NPP 表现 为 显著 负 相 关 , 且 存在 3 个 月 的 滞后 效 
应 。 暖 夜 日 数 TN90P 与 春季 同期 NPP 表 现 为 极 显 
著 正 相关 ,存在 2 个 月 的 滞后 效应 。 春 季 平 均 温 度 
TEM 与 同期 NPP 表 现 为 极 显 著 正 相关 , 且 存 在 3 个 
月 的 滞后 效应 。 总 体 而 言 ,春季 均 温 升 高 以 及 暖 夜 
日 数 增加 均 有 利于 植被 NPP 增 加 。 而 夏季 DTR HE 


在 裸露 区 ,NPP 主要 受 GSL , PRCPTOT 和 日 最 低 气 
温 极 大 值 TNx 的 影响 。 覆 沙 区 主要 受到 DTR RE 
日 数 TX90P、TNx 的 影响 。 履 土 区 主要 受到 GSL, 
DTR 和 TX90P 的 影响 。 

2.6 气候 和 人 类 活动 对 而 砂岩 区 年 NPP 变 化 的 相 


结合 极端 气候 对 NPP 影 响 的 重要 性 排序 结 
(图 8) ,得 到 气温 日 较 差 DTR 对 全 域 NPP 变 化 的 重 
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* 表 示 在 0.05 水平 ( 双 侧 ) 显 著 ;** 表 示 在 0.01 水 平 ( 双 侧 ) 显 著 。 
滞后 1 个 月 .2 个 月 和 3 个 月 NPP 的 相关 性 


Fig. 7 Correlation coefficients of precipitation and temperature index with NPP for time lags in the Pisha sandstone area 
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Fig.8 Importance of 18 extreme climate indices to NPP in various regions of the Pisha sandstone area 


要 性 最 高 ,因此 ,本 研究 利用 多 元 线性 回归 构建 NPP 
与 年 均 温 度 .年 降水 量 以 及 气温 日 较 差 的 关系 , 进 
行 多 元 线性 回归 残 差 和 相对 贡献 分 析 , EEDA 
区 NPP 变 化 的 驱动 因素 情况 (图 9) ,结果 显示 ,气候 
因素 和 人 类 活动 共同 作用 是 而 砂岩 区 年 NPP 增 加 
的 主要 原因 ,面积 约 占 93.25%。 和 气候 因素 和 人 类 活 


促进 作用 。 

进一步 分 析 气 候 因 素 和 人 类 活动 对 而 砂岩 
NPP 增 加 的 相对 贡献 (图 10) ,结果 显示 ， 
全 域 平均 气候 贡献 度 为 50.23% ,平均 人 类 活动 贡献 
率 为 49.77% ,气候 和 人 类 活动 对 NPP 增 加 贡献 相 
365 (ACE RIAD AAS [8] DX Sak TE SF LER UE AS Fn] ,裸露 


动 共 同 导 致 NPP 减 少 的 面积 约 占 1.7% , 仅 受气 候 因 
素 积极 影响 的 区 域 占 比 为 2.7% , FALL, aL eb A AL 
被 NPP 变化 主要 受到 气候 因素 和 人 类 活动 的 积极 


区 和 和 窗 沙 区 气候 贡献 占 主导 地 位 (相对 贡献 率 超过 
5096) ,覆土 区 人 类 活动 占 主导 地 人 位。 裸露 区 平均 气 
候 贡 献 率 为 62.13%, 人 类 活动 贡献 率 为 37.87%。 禾 
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Fig.9 Spatial distribution of NPP drivers in vegetation in the 
Pisha sandstone area from 2001 to 2021 


沙 区 平均 气候 贡献 率 为 60.06%, 人 类 活动 贡献 率 为 
39.94%。 履 土 区 平均 气候 贡献 率 为 39.60%, 人 类 活 
BH vt AA y 60.40% . 3X 5j Vr AE A zs DX — RY 
生态 恢复 措施 的 实施 密切 相关 ,通过 治理 水 土 流 
失 , 而 砂岩 区 植被 覆盖 度 增加 ,植被 生长 状态 改 
善 。 旦 相 较 于 裸露 区 和 覆 沙 区 ,覆土 区 人 类 活动 对 
NPP 增 长 贡献 起 主导 作用 ,植被 恢复 工程 实施 效果 
更 为 显著 。 


3 讨论 


3.1 Aiba tee NPP 时空 变化 

2001—2021 年 而 砂岩 区 NPP 整 体 呈 波动 上 升 
趋势 ,其 中 ,和 窗 土 区 NPP 增 长 速率 高 于 窗 沙 区 和 裸 
露 区 ,这 是 由 于 覆土 区 土壤 养分 含量 高 "", 具 有 较 
好 的 立地 条 件 , 生 态 恢复 工程 效益 在 覆土 区 效益 显 
著 , 人 类 活动 对 覆土 区 NPP 变化 的 贡献 达 60.4%， 
此 ,NPP 增 速 最 快 。 而 受 立地 条 件 等 因素 的 制约 , 覆 
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气候 贡献 率 /% 
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沙 区 和 裸露 区 植被 存活 率 较 低 ” ,植被 建设 达到 的 
效果 相对 有 限 汪 ,人 类 活动 对 NPP 增 加 的 贡献 小 于 
气候 因素 。 从 全 域 整 体 来 看 ,人 类 活动 和 气候 因素 
对 NPP 增 长 的 贡献 相当 ,说 明 而 砂岩 区 生态 恢复 措 
施 实 施 效果 显著 ,但 未 来 NPP 由 增加 变 为 减少 的 区 
域 达 82.50%, 植 被 退化 风险 较 大 ,这 与 Chen 等 对 
黄土 高 原 植被 研究 结果 一 致 , 即 近 21 a 来 黄土 高 原 
NPP 有 所 增加 ,但 在 未 来 很 可 能 无 法 持续 增加 。 
此 ,还 需要 继续 实施 生态 环境 保护 工程 ,调整 未 来 
植被 恢复 策略 ,巩固 大 t 砂 岩 区 植被 生态 恢复 效益 ， 
促进 植被 恢复 建设 可 持续 发 展 。 
3.2. SIRARI IA NPP 的 影响 

年 际 尺度 上 ,极端 冷 .极端 暖和 气温 日 较 差 
DTR 25] 5j fit A [X NPP 为 负 相 关 , 这 与 刘铮 等 以 
黄土 高 原 为 研究 区 的 研究 结果 一 致 。 而 年 均 温 度 
TEM 生长 季 长 度 GSL 与 NPP 表 现 为 正 相 关 关 系 ,说 
明年 均 温 增加 ,可 能 通过 增加 植被 生长 季 长 度 来 促 
进 植被 生长 发 育 ”” ,但 温度 存在 阅 值 ”, 当 温度 过 
高 时 ,导致 植被 呼吸 作用 增强 , 植物 气孔 关闭 进而 
影响 植被 光合 作用 ,温度 过 低 时 ,植物 生长 所 需 
的 积温 不 足 ,植被 发 育 缓慢 ,不 利于 植被 生长 ”。 
气温 日 较 差 DTR 5 lb A DX NPP 主要 表现 为 显著 
负 相 关 , 这 与 何 航 等 "研究 发 现 ,中 国 北方 生长 季 
NDVI 与 DTR 呈 负 相 关 的 结论 一 致 ;同时 ,Xu 等 外 也 
人 研究 发 现 , 中 国 北方 沿海 地 区 NPP 与 DTR 5 fA 4H 
关 。 可 能 是 一 方面 ,DTR 受 温 度 波动 大 小 的 影响 ， 
会 对 植被 强制 休眠 产生 作用 ,影响 发 芽 的 频率 1。 
另 一 方面 , Wan 等 “通过 在 中 国 北方 温带 草原 增 温 
试验 研究 发 现 , 白 天 增 温 不 利于 植被 生长 ,夜间 增 
温 可 以 促进 植被 生产 力 提 高 ,DTR 是 日 最 高 气温 和 
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图 10 tee LZR ALA AST SIDE LAD t DX NPP 的 相对 贡献 


Fig. 10 Relative contributions of climatic factors and human activities to NPP in the Pisha sandstone area 
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日 最 低 气 温 的 差 值 , 受 极端 气温 的 影响 ,DTR 减 小 
即日 最 高 气温 降低 或 日 最 低 气 温 升 高 。 因 此 ,DTR 
减 小 可 能 有 助 于 而 砂岩 区 植被 生长 。 相 比 之 下 slt 
砂岩 区 年 NPP 与 极端 强 降 水 指数 均 整 体 表 现 为 正 
相关 关系 ,这 与 任 晋 媛 等 ”人 研究 发 现 ,内 蒙古 NPP 
与 极端 降水 均 呈 正 相关 的 结论 一 致 , 这 是 由 于 而 砂 
宕 区 属于 半 干 旱地 区 ,植被 生长 受降 水 量 限 制 , 且 
植被 恢复 导致 区 域 径流 和 土壤 可 用 水 量 进一步 减 
^p ges hk s PRG, H BID ALL Ab S Dx A FE 
水 量 增多 和 降水 强度 增 大 有 利于 植被 NPP 增 长 。 
季 方 尺度 上 ,春季 均 温 以 及 冷 极 值 温度 升 高 均 有 助 
于 植被 生长 ,这 可 能 与 春季 增 温 促进 植被 生长 季 提 
前 有 关 ” ,而 夏季 高 温 与 干旱 不 利于 NPP 增 加 , 干 
旱 会 影响 植被 光合 作用 和 呼吸 作用 ,导致 植被 生长 
能 力 下 降 ,削弱 碳 汇 功 能 ”。 
3.3 人 类 活动 对 础 砂岩 区 NPP 的 影响 

人 类 活动 通过 改变 础 砂岩 区 土地 利用 类 型 ”， 
改良 砚 砂 岩 区 土壤 质地 ,改善 立地 条 件 ” ,影响 林 
草 等 植被 生长 ,促进 区 域 植被 恢复 。 目 前 aI E 
区 以 沙棘 生物 措施 为 治理 突破 口 ,实施 了 一 系列 的 
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生态 环境 治理 工程 。1998 年 水 利 部 正式 实施 晋 陕 
ex IER ey CUP RAE AS TFE , 2008 EF JE [Bee lE 
A DX fe IF TR Di o Ub LE AS AD TE , 2013 年 十 大 孔 
兑 沙 环 生 态 减 沙 工 程 实施 ”。 而 砂岩 区 已 引入 大 
量 水 土 保 持 植 被 如 油 松 PR PARE SE T E 
岩 区 植被 覆盖 度 光 ,改善 了 水 土 流失 的 程度 3。 
2001—2021 年 间 而 砂岩 区 林地 和 建设 用 地 面积 增 
加 ,草地 耕地 和 未 利用 土地 面积 减少 ( 表 3) ,林地 
MEARE A TH, MAKERE 
土地 利用 类 型 转换 的 面积 占 总 面积 的 25.12% ,土地 
利用 变化 区 域 促 进 NPP 增 加 量 仅 占 总 NPP 增 加 量 
的 24.58%( 表 4)。 说 明 而 砂岩 区 NPP 增 加 除了 与 十 
地 利用 类 型 转换 有 关 ,更 与 未 发 生 土 地 利用 变化 区 
域 的 植被 生长 状况 改善 有 关 。 近 20 a 来 , 随 着 退耕 
还 林 、 退 牧 还 草 等 政策 的 实施 , 而 砂岩 区 退化 的 草 
地 得 到 了 和 较 好 的 恢复 ,草地 NPP 增 加 。 且 人 工种 
植 的 沙 杯 和 油 松 对 而 砂岩 区 土壤 碳 、 氮 含量 的 提升 
BR BENE , 促 进 了 而 砂 岩 区 林地 NPP 增 加 。 人 类 
活动 对 楷 砂岩 区 植被 NPP 以 促进 作用 为 主 ,这 与 已 
Ae S077 REA AL a KR STE 


#3 而 砂岩 区 2001 一 2021 年 土地 利用 转移 矩阵 


Tab.3 Transition matrix of land use types in the Pisha sandstone area between 2001 and 2021 I% 
2001 年 - EMT 转 出 总 量 
草地 建设 用 地 耕地 林地 水 域 未 利用 土地 
草地 52.94 4.36 4.67 1.47 0.36 0.93 11.80 
建设 用 地 0.19 0.97 0.09 0.02 0.02 0.01 0.33 
耕地 5.39 1.70 13.19 0.47 0.26 0.11 7.93 
林地 0.70 0.36 0.22 2.63 0.04 0.03 1.35 
水 域 0.48 0.23 0.30 0.09 1.86 0.04 1.14 
未 利用 土地 2.02 0.27 0.12 0.13 0.04 3.28 2.57 
转 人 总 量 8.78 6.92 5.40 2.18 0.72 1.12 25.12 
#4 WARK 2001—2021 ib Fl AZE NPP Æ rh E 
Tab. 4 Proportion of NPP transferred during the change in land use types in the 
Pisha sandstone area between 2001 and 2021 /% 
2001 年 EUR 转 出 总 量 
草地 建设 用 地 耕地 林地 水 域 未 利用 土地 
草地 55.48 2.91 5.48 1:22 0.33 0.91 10.85 
建设 用 地 0.19 0.81 0.09 0.02 0.01 0.01 0.32 
耕地 6.41 1.05 14.70 0.39 0.25 0.11 8.21 
林地 0.93 0.30 0.30 2.64 -0.05 0.04 1.52 
水 域 0.45 0.20 0.27 0.07 1.74 0.03 1.02 
未 利用 土地 2.23 0.19 0.10 0.10 0.04 0.05 2.66 
转 人 总 量 10.21 4.65 6.24 1.8 0.58 INI 24.58 
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对 区 域 生 态 环境 建设 发 挥 了 重要 作用 。 


4 结论 

本 文 基于 CASA BRAY TE AID A XA BRI 
生产 力 , 分 析 其 多 年 时 空 变化 特征 ,结合 年 均 温 度 
和 18 个 极端 气候 指数 分 析 而 砂 岩 区 NPP 对 气候 变 
化 的 响应 ,并 进一步 计算 得 到 气候 和 人 类 活动 对 
NPP 的 相对 贡献 。 具 体 结论 如 下 : 

(1) 空间 分 布 上 , RIERP KX NPP 大 小 与 覆土 程 
度 分 布 一 致 , 即 NPP 值 在 覆土 区 AWK .裸露 区 逐 
HEIR 2001—2021 年 间 , MEEI KER NPP E fe 
著 增加 趋势 ,但 各 区 域 NPP 增 长 均 没 有 较 强 持续 
PE, TEARS ALE KAR 82.59 HKI NPP FH 4 in 
RS Ag li > SN, BE Db est DX AS HA AS TG IB HY) 
风险 。 

(2) 而 砂 岩 区 NPP 对 各 气候 指标 的 响应 不 同 。 
年 际 尺 度 上 ,与 NPP 相 关 性 强 的 气温 指标 主要 是 气 
温 日 较 差 DTR .生长 季 长 度 GSL EA TEM 27 
夜 日 数 TNI0P。NPP 与 TEM GSL 表现 为 广泛 的 正 
相关 ,而 与 TN10P、DTR 表现 为 负 相 关 。 相 反 ,不 同 
强度 的 极端 降水 指数 (PRCPTOT、RX1day , RX5day , 
R95P、R99P、SDIT .CWD) 与 年 NPP 均 为 正 相 关 关 
系 ,表明 当前 极端 强 降水 对 夏 砂岩 区 NPP 起 积极 影 
响 。 季 节 尺 度 上 ,春季 和 夏季 NPP 对 平均 气候 态 与 
极端 气候 的 响应 具有 明显 的 滞后 效应 ,春季 增 温 有 
助 于 NPP 增 加 ,夏季 降水 增多 显著 促进 植被 生长 ， 
但 夏季 暖 凡 日 数 增 加 不 利于 植被 生长 。 

(3) 近年 来 而 砂 岩 区 一 系列 生态 恢复 工程 的 实 
施 显著 促进 了 植被 生长 ,人 类 活动 和 气候 变化 对 全 
域 NPP 增 长 的 贡献 相当 ,分 别 为 49.77% 和 50.23% 。 
各 区 域 上 略 有 不 同 ,裸露 区 和 禾 沙 区 NPP 变化 以 气 
候 因 素 为 主导 ,气候 因素 贡献 率 分 别 为 62.13% 和 
60.06%。 覆 土 区 NPP 变化 以 人 类 活动 为 主导 ,人 类 
活动 页 献 率 为 60.4%。 

(4) 础 砂岩 区 植被 生长 受气 候 因 素 和 人 类 活动 
的 共同 影响 ,促进 而 砂岩 区 植被 稳定 健康 生长 ,应 
持续 推进 林 草 植被 保护 与 建设 工作 ,加强 土壤 侵蚀 
治理 ,优化 覆土 区 植被 格局 ,充分 发 挥 覆 土 区 植被 
恢复 潜力 ,加 强 对 裸露 区 和 履 沙 区 植被 保护 工作 ， 
并 通过 人 为 措施 缓解 夏季 干旱 与 高 温 对 植被 生长 
的 消极 影响 。 
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Impact of climatic factors and human activities on the net primary productivity 
of the vegetation in the Pisha sandstone area 


WANG Yi'en", RAO Liangyi"" 
(1. College of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China; 2. Key Laboratory of State 
Forestry and Grassland Administration on Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China) 


Abstract: As a highly serious soil and water loss area in the Loess Plateau, vegetation growth plays an important 
role in controlling soil erosion and maintaining ecological balance. This study calculated the net primary 
productivity of vegetation NPP and 18 extreme climate indices in the Pisha sandstone area from 2001 to 2021 
based on the CASA model and Rclimdex 1.0, respectively. Trend analysis, correlation analysis, random forest 
importance ranking, and residual analysis were used to study the spatial- temporal variation of NPP and its 
response to climate factors in the Pisha sandstone area. The relative contributions of climate factors and human 
activities to the NPP in the Pisha sandstone area were also calculated. The results showed that (1) NPP variation 
in all regions of the Pisha sandstone area from 2001 to 2021 had a significant increasing trend, but in the future, 
82.5% of the NPP in the Pisha sandstone area will change to a decreasing trend. (2) On the annual scale, NPP 
correlated positively with average annual temperature, total annual precipitation, and extreme heavy precipitation 
index and correlated negatively with cold night days TN10P and diurnal temperature range DTR. On the seasonal 
scale, the increase in average temperature and warm night days in spring was conducive to increase NPP, and 
there is a lag effect. Increasing the number of warm days in summer was unconducive to vegetation growth, and 
the NPP has a three-month lag response to the number of warm days in summer. Extreme heavy precipitation in 
summer was conducive to NPP increase, whereas summer drought was unconducive to vegetation growth, and 
NPP has a three-month lag response to the number of continuous dry days. (3) Both climate change and human 
activities contribute positively to NPP in the Pisha sandstone area. The climate contribution of the bare area and 
the covered sand area is dominant (62.13% and 60.06%, respectively), whereas that of the covered soil area is 
dominated by human activities (60.40%). 

Keywords: Pisha sandstone area; net primary productivity of vegetation(NPP); extreme climate; human 


activities; residual analysis 


